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Abstract. Magnetic field-aligned electric fields are character- 
istic features of magnetic reconnection processes operating in 
externally agitated magnetized plasmas. An especially interest- 
ing environment for such a process are the coronae of accretion 
disks in active galactic nuclei (AGN). There, Keplerian shear 
flows perturb the quite strong disk magnetic field leading to in- 
tense current sheets. It was previously shown that given field 
strengths of 200 G in a shear flow, reconnection driven mag- 
netic field aligned electric fields can accelerate electrons up to 
Lorentz factors of about 2000 in those objects thus providing 
us with a possible solution of the injection (pre-acceleration) 
problem. However, whereas in the framework of magnetohy- 
drodynamics the formation of the field-aligned electric fields 
can be described consistently, the question has to be addressed 
whether the charged particles can really be accelerated up to the 
maximum energy supplied by the field-aligned electric poten- 
tials, since the accelerated particles undergo energy losses ei- 
ther by synchrotron or inverse Compton mechanisms. We pre- 
sent relativistic particle simulations starting from electric and 
magnetic fields obtained from magnetohydrodynamic simula- 
tions of magnetic reconnection in an idealized AGN config- 
uration including nonthermal radiative losses. The numerical 
results prove that the relativistic electrons can be effectively 
accelerated even in the presence of an intense radiation bath. 
Energies from 50MeV up to 40GeV can be reached easily, de- 
pending on the energy density of the photon bath. The strong 
acceleration of the electrons mainly along the magnetic field 
lines leads to a very anisotropic velocity distribution in phase 
space. Not even an extremely high photon energy density is 
able to completely smooth the anisotropic pitch angle distri- 
bution which is characteristic for quasi monoenergetic particle 
beams. 
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1. Introduction 

Active galactic nuclei (AGN) can be regarded as accreting su- 
permassive black holes surrounded by accretion disks (Camen- 
zind 1990; Miyoshi et al. 1995; Burke and Graham-Smith 1997 
and references therein). Relativistic electrons in AGN reveal 
themselves by hard X-ray (probably due to pair production (cf. 
Svensson 1987; Done and Fabian 1989)) and 7-ray emissions 
as well as radio observations of superluminous motions (e.g. 
Abraham et al. 1994). The 7-radiation observed from quasars 
and blazars may originate in a distance R of 10 2-3 gravita- 
tional radii from the central engine (Dermer and Schlickeiser 
1993). At that distance no significant pair production happens 
but relativistic leptons scatter via the inverse Compton process 
the IR-UV radiation of the accretion disk within a relativisti- 
cally moving jet. It is well known that "standard" mechanisms 
for the acceleration of high energy leptons, as diffusive shock 
wave acceleration and resonant acceleration by magnetohy- 
drodynamic al (MHD) turbulence can only work efficiently for 
Lorentz factors larger than 7 cr i t ~ m p /m c (where m p and m e 
denote the proton and electron masses). Consequently, charged 
particles accelerated via shocks or MHD turbulence have to be 
pre-accelerated which confronts us with the injection problem 
(Blandford 1994; Melrose 1994) in the AGN context. 

In a differentially rotating magnetized accretion disk gas, 
the evolution of a magnetized corona is quite hard to sup- 
press (Galeev et al. 1979; Stella & Rosner 1984). Driven by 
the buoyancy force magnetic flux tubes ascend into the disk 
corona, thereby their footpoints are sheared by the differential 
rotation of the disk. Either by internal shear or by encounter- 
ing already present magnetic flux, magnetic reconnection and 
accompanied rapid dissipation of magnetic energy happens in 
the coronal plasma. Such a behavior can be studied with great 
detail for example in the solar corona (Parker 1994 and refer- 
ences therein). In a recent contribution we investigated the pos- 
sible role of magnetic field-aligned electric fields (En) in the 
context of magnetic reconnection operating in AGN coronae 
for the pre-acceleration of leptons (Lesch and Birk 1997; here- 
after LB). It could be shown that field-aligned electric potential 
structures in relatively thin current sheets form. For reasonable 
physical parameters such electric potentials are strong enough 
to accelerate electrons up to 7 w 2000, in principle. However, 
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in the framework of MHD the actual energies of the acceler- 
ated particles cannot be calculated. It is the aim of the present 
contribution to corroborate the model introduced in LB with 
the help of relativistic particle simulations by taking macro- 
scopic electric and magnetic field configurations obtained by 
the MHD simulations as an input. 

In the next section we resume the MHD model in a nutshell 
and present the details of the resulting three-dimensional elec- 
tric and magnetic fields. In Sec. |]we discuss our approach to 
the numerical study of high-energy particles and show the nu- 
merical results dwelling on particle spectra and energies. Even- 
tually, we discuss our findings in Sec. Q 



2. The MHD model for magnetic reconnection in the 
corona of an accretion disk 

The MHD framework for the particle simulations is modeled 
by means of a 3D resistive MHD code that integrates the one- 
fluid balance equations (e.g. Krall and Trivelpiece 1973): 
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where p, v, p, B, j and 77 denote the plasma mass density, veloc- 
ity and pressure, the magnetic field, the electric current density 
and the electrical resistivity. 

The scenario we have in mind can be briefly described as 
follows (for details see LB): Different convective plasma mo- 
tions at different regions of an AGN corona or differential shear 
motion of the disk result in a sheared magnetic field and conse- 
quently, in the formation of field-aligned electric current sheets. 
When this convective shear motion is strong enough to result 
in a supercritical current density (drift velocity larger than the 
local sound velocity c s = ^/k^,T c /m v current-driven microin- 
stabilities are excited which lead to anomalous electrical resis- 
tivity (e.g. Huba 1985; Benz 1993). The formation of localized 
regions of anomalous dissipation or, more general, the local 
violation of ideal Ohm's law, can be regarded as the neces- 
sary onset condition for magnetic reconnection (Schindler et 
al. 1988, 1991). During the dynamical evolution of the recon- 
nection process very localized acceleration regions (E|| ^ 0) 
form. 

Since we want to concentrate on a single sheared coronal 
loop which is represented by a current-carrying magnetic flux 
tube we start from an appropriate idealized initial configuration 
for the MHD simulation run which is characterized by a homo- 




Fig. 1. The x-,y- and z-component (from top to bottom) of the 
reconnection magnetic field at the height of the central recon- 
nection region after t = 120ta. 



geneous plasma and a force-free magnetic field (cf. Birk and 
Otto 1996; LB): 



B = B y0 tanh(a;) e v - * / B z0 2 
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where B z0 = 5 and B y0 = 1 denote the normalized con- 
stant main component and the shear (toroidal) component of 
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Fig. 2. The y- and z -component of the reconnection magnetic 
field at the y = O-half plane (i.e. the middle on the of the central 
reconnection region) after t = 120ta- 
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the magnetic field, respectively. As discussed in LB we choose 
for the half-width of the current sheet w as 10 7 cm, a shear 
magnetic field of B y o = 200G, a particle density of n = 
3 • 10 5 cm -3 and 10 4 statamp cm -2 for the critical current den- 
sity of jcrit, which corresponds to an electron temperature of 
6-10 r K (cf.Ulrich 1991;Nandra and Pounds 1994). The mag- 
netic Reynolds number is chosen as S = AirvAw/r/c 2 = 10 3 . 

The initial equilibrium configuration is perturbed by some 
sheared flow: 



0.02 



Vy(X,Z,t) 
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with an amplitude chosen as v y o = 0.5% of the shear Alfven 
velocity va — B y o/ This perturbation is transported 
down along the main component of the magnetic field via shear 
Alfven waves during the dynamical evolution. An anomalous 
resistivity will be switched on when the current density exceeds 
the critical value j C rit- The resistivity is localized in height 
along the poloidal component of the magnetic field in order 
to model an acceleration region of the length of 40iu as well as 
in the y-direction (for technical details we again refer to LB). 
During the dynamics reconnection electric and magnetic fields 




Fig. 3. The x-,y- and z-component (from top to bottom) of the 
reconnection electric field after t — 120ta- 



evolve which can be used as an input configuration for the par- 
ticle simulation of the acceleration of high-energy electrons in 
the reconnection region. A snapshot after t = 120 dynamic 
times (i.e. ta = w/va) of the Cartesian components of the 
magnetic and electric fields are shown in Figures The re- 
sulting x-component of the magnetic field shows the typical 
bipolar reconnection characteristic (Fig. |IJ top panel). The mid- 
dle panel of Fig. |l] shows the y-component of the magnetic field 
associated with the current sheet and the lower panel shows the 
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U, t irne 



= 120 




Fig. 4. The generalized electric potential U after t — 120ta- 



z-(poloidal) component. Fig. |] shows the y- and z-component 
of the magnetic field at the y = O-half plane. The y-component 
is enhanced locally due to the sheared flow applied at the upper 
boundary and transported along the poloidal component of the 
magnetic field via shear Alfven waves. The three-dimensional 
electric vector field at the height of the central reconnection re- 
gion is illustrated in Fig. |[ The convective x-component (top 
panel) is the result of the applied sheared flow and the diver- 
gent reconnection flow whereas the convective y-component 
(middle panel) is due to the convergent (v x (x)) reconnection 
flow. The dominant accelerating component is E z (lower panel) 
which is caused by the finite anomalous electrical resistivity 
E z ~ r)j z and accordingly localized in the x- and y-direction. 

The generalized field-aligned electric potential in normal- 
ized units U = — J Ends (ds is a magnetic field line element) 
associated with the electric and magnetic fields is illustrated in 
Fig. Q It is this elongated relatively thin potential structure by 
which the electrons are accelerated as will be discussed in the 
next section. 



3. Particle simulations 

For the sake of simplicity, we perform our simulations as rel- 
ativistic test particle simulations. This means, that any inter- 
action between the simulated particles are neglected, as well 
as the backreaction of the considered particles on the electro- 
magnetic fields. In our calculations we use the stationary fields, 
which are explained in Sec. ||, simply as the "stage" for the ac- 
celerated charged particles. 

In addition to the acceleration force, we consider the radia- 
tive losses in our simulations, especially the synchrotron (SY) 
and the inverse Compton (IC)-losses. In particular, we are in- 
terested in the question: 

Is magnetic reconnection able to accelerate particles to rel- 
ativistic energies in the presence of synchrotron and inverse 
Compton losses? 



3.1. Momentum Equation and Numerical Procedure 

From the mathematical point of view the problem of test parti- 
cle simulation is an integration of the equations of motion for a 
single charged particle. Since we are dealing with test particles, 
we have to integrate: 

dr 

dt 
dp 

dt 

This simplifies our problem enormously. As usual, the vec- 
tors r, v and p represent position, velocity and momentum of 
the particle and F is the force acting on the particle. We use the 
relativistic relations 
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between velocity and momentum, with the rest mass of the par- 
ticle too and the velocity of light c. 

The resulting force acting on the particle is, in the simplest 
case, given by the Lorentz force Fl = q [E + v/c x B] , where 
q denotes the electric charge of the particle. This is used in cal- 
culations in which radiative loss processes are neglected. If we 
take these losses into account, we require in addition a radia- 
tive loss force, which describes the corresponding interaction 
between radiation field (photons) and particles. Such a force 
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can be found in Landau and Lifshitz (1951). It consists of every 
radiative loss process due to field-particle interaction, includ- 
ing SY- and IC-losses. Losses due to particle-plasma-wave 
interaction are not considered up to now. 

The force given in Landau and Lifshitz (1951) is used here 
in a modified form. All terms with temporal or spatial deriva- 
tives are neglected, because the external fields are varying on 
great scales and it is very expensive to calculate those terms for 
the microscopic photon bath. So we consider the IC-losses by 
including an additional term, which describes the space- and 
time-averaged influence of the IC-losses of an isotropic pho- 
ton bath with the energy density [/R a d- After some calculations 
this leads to the modified expression 
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In this equation (3 = v/c is the velocity normalized to the 
speed of light. 

The choice of numerical routines needs some comments. 
The routine, which is executed most, is the evaluation of the 
force F = Fl + FR a( j at every timestep. This evaluation con- 
sists of an interpolation of the fields, which are known only on 
a grid. For the sake of rapidness we use a linear interpolation. 
This means, that at the grid boundaries the derivatives of our 
field quantities are not defined. For this reason we use a Runge- 
Kutta procedure of second order accuracy in combination with 
an adaptive stepsize control. 

3.2. Particle simulations: Energies and Spectra 

In our calculations we have obtained the following results. 
The SY-losses are completely unimportant in our AGN con- 
text. In Fig. |^ one can see clearly, that the curves without any 
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Fig. 5. Temporal evolution of the kinetic energy of a specific 
electron with respect to different loss processes and strengths 

losses and with SY-losses only are nearly identical. In the fol- 
lowing all energies are given in units of the rest energy of 
the electrons m c c 2 and [/Mag denotes the energy density of 
the magnetic field. For IC-losses the situation is completely 
different. The final energies decrease steadily with increasing 
photon energy density. However, even at the extreme value of 
C/R a d — 100 • C/Mag the electrons gain some several tens of 
MeV. Anyway, one can clearly see that for reasonable pho- 
ton energy densities ([/Mag — [^Rad) a significant part of the 
injected electrons are accelerated towards GeV-energies. This 
can solve the injection problem. 

The final energy spectra of the accelerated electrons for dif- 
ferent strengths of IC-losses are shown in Fig. ||. For all spec- 
tra, the electrons were injected homogeneously distributed in 
the acceleration region with a velocity-distribution that cor- 
responds to a Maxwell-distribution of temperature 6 ■ 10 8 K. 
In any case, the resulting final spectrum is a very broad one, 
which goes over some decades. The effects of IC-losses are 
clearly visible. The corresponding maximum of the distribu- 
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Fig. 6. Final energy distributions of the accelerated electrons 
for three different IC-loss strengths (top: no losses, mid- 
dle: Ubm = 10 ■ [/Mag, bottom: [/ Rad = 100 ■ U Mag ) 



tions is shifted towards lower energies with increasing photon 
densities. It shifts from about 8GcV without energy losses to 
approximately lGeV at an energy density of lOOC/Mag- The 
same tendency holds for the maximum energy achievable by 
the electrons. The maximum energy gain without losses coin- 
cides with the energy gain one would expect from the potential 
shown in Fig. Due to the fact, that the minimum energy gain 
in these distributions is quite constant with respect to a vari- 
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Fig. 7. Comparison between the distribution of the final mo- 
mentum parallel (top) and perpendicular (bottom) to the mag- 
netic fieldlines (Note that the abscissae are different in the two 
figures) 

ation of E/r< a d (this results from the 7 2 -dependence of FR a( j), 
the distribution becomes steeper and narrower with stronger en- 
ergy losses. 

Fig. [7] shows the distributions of the final momenta of the 
particles perpendicular and parallel to the magnetic field. It is 
evident, that the motion of the electrons is strongly anisotropic. 
They are mainly following the magnetic field lines, whereas 
any motion vertical to the field lines is suppressed. Thus, one 
may conclude that the accelerated electrons represent a collec- 
tion of quasi-monoenergetic beams. It is not shown, that the 
movement becomes the more isotropic, the higher the losses 
are. But this is straightforward, because in the loss force used 
FRad we have assumed an isotropic photon bath which natu- 
rally tends to isotropize the electron orbits. 

In Fig. ^ and Fig. [7] (top) one also can see another char- 
acteristic feature of our field configuration. It is quite obvious 
that in all four distributions there is something like an absorp- 
tion line between the two local maxima. This effect can be ex- 
plained in a natural way, if one looks directly at some chosen 
trajectories of the particles. Fig. [| shows the trajectories of the 
electrons in combination with a projection of the orbits on the 



Fig. 8. Trajectories of a chosen set of electrons. The projections 
on the walls behind, were added artificially for a better 3-D 
impression. 

walls behind. These projections show, that there are two popu- 
lations of electrons. One population is accelerated towards and 
into the central region of reconnection, whereas the other pop- 
ulation is repelled. A comparison of the final energies of the 
electrons with the membership in one of the two classes shows 
that exactly these electrons which are forced into the reconnec- 
tion zone are accelerated the most. The members of the other 
group of electrons only gain some moderate amount of energy. 
This indicates that the instability of the orbits is the cause for 
the appearance of the "absorption line". This effect is quite in- 
teresting, because it represents a selection process, which leads 
to an even stronger beaming of the attracted electrons and to an 
even stronger anisotropy. 

4. Summary and Discussion 

We addressed the question of charged particle acceleration dur- 
ing magnetic reconnection processes by means of relativistic 
test particle simulations. Whereas this point is crucial for a 
great variety of cosmic plasmas in this contribution we dwelled 
on the pre-acceleration problem in the AGN context. Possible 
further applications include among others as different plasma 
systems as the terrestrial discrete auroral arcs, the solar coronal 
flares, radio activity in T-Tauri magnetospheres, non-thermal 
emission at the edges of high-velocity clouds that hit the galac- 
tic plane and the generation of electron beams in the magneto- 
spheres of neutron stars. The starting point for our present in- 
vestigations are results obtained by an MHD simulation study 
(LB). In contrast to previous work carried out by different 
other groups, we were able to study particle acceleration in 
large-scale non-linearly developed reconnection electromag- 
netic fields rather than being restricted to the prescription of 
somewhat idealized analytical field solutions. Moreover, since 
for the considered parameter regime we have to expect an 
intense radiation field the test particle simulations were per- 
formed including the relevant radiative losses. A very impor- 
tant question is whether charged particles can really be accel- 



S chopper et al. 



7 



erated in the reconnection region up to the pretty high energy 
values one might deduce from the fluid treatment. Our find- 
ings indicate that, in fact, as expected from the fluid simulations 
(LB), leptons can be accelerated in reconnection zones located 
in AGN coronae up to high Lorentz factors; i.e. a significant 
portion of test particles gain about the maximum energy pro- 
vided by the generalized field-aligned electric potential struc- 
tures formed during the magnetic reconnection processes we 
have modeled within the framework of MHD. Thus, particle 
acceleration in reconnection zones may be considered as a way 
out of the injection problem we face in the AGN context. 

Since sheared magnetic fields can be expected as a very 
common phenomena in cosmic environments like accretion 
disks, stellar coronae and interstellar medium, we think that 
our relativistic particle studies can be regarded as realistic with 
respect to the used configuration and the dominant forces. In 
this contribution the magnetic reconnection process is driven 
by some resistive mechanism originating from plasma instabil- 
ities. The excited electromagnetic oscillations serve as resis- 
tance. We note that in plasmas which are collisionless (both no 
Coulomb collisions and no turbulent wave excitation), the par- 
ticle inertia presents the ultimate source of resistivity and for 
the magnetic dissipation. The sheared magnetic fields in colli- 
sionless systems evolve into very thin filaments, in which the 
lifetime of the particle determines the electrical conductivity, 
thereby allowing for efficient dissipation via effective particle 
acceleration (Lesch and Birk 1998). 

Our simulations are test-particle simulations, thus, we 
plan for future studies to include, additionally, ponderomotive 
forces and the back reaction of the current carried by the high 
energy particles. Whether or not the latter aspect becomes im- 
portant depends on the density of the run-away electrons lim- 
ited by the Dreicer electric field (e.g. Benz 1993). 
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